Noise of a single-electron transistor in the regime of large quantum fluctuations of 

island charge out of equilibrium 



O 
O 

> 
O 



Yasuhiro Utsumi,* Hiroshi Imamura, Masahiko Hayashi, and Hiromichi Ebisawa 
Graduate School of Information Sciences, Tohoku University, Sendat 980-8579, Japan 

(Dated: February 1, 2008) 

By using the drone-fermion representation and the Schwinger-Keldysh approach, we calculate the 
current noise and the charge noise for a single-electron transistor in the non-equilibrium state in 
the presence of large quantum fluctuation of island charge. Our result interpolates between those 
of the "orthodox" theory and the "co-tunneling theory". We find the following effects which are 
not treated by previous theories: (i) At zero temperature T = and at finite applied bias voltage 
eV^j ^ Tk, where Tk is the "Kondo temperature", we find the Fano factor is suppressed more than 
the suppression caused by Coulomb correlation both in the Coulomb blockade regime and in the 
sequential tunneling regime, (ii) For T ^ \eV\/2 ^ Tk, the current noise in the presence of large 
charge fiuctuation is modified and deviates from the prediction of the orthodox theory. However, 
the Fano factor coincides with that of the orthodox theory and is proportional to the temperature, 
(iii) For eV, T < Tk, the charge noise is suppressed due to the renormalization of system parameters 
caused by quantum fluctuation of charge. We interpret it in terms of the modification of the "unit" 
for island charge. 

PACS numbers: 73.23.Hk,72.70.-|-m 
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I. INTRODUCTION 



In a small metallic island where the charging en- 
ergy Ec exceeds the temperature T (we use the unit 
/cb = 1), the Coulomb interaction affects transport prop- 
erties through the island. The resulting phenomenon is 
called the Coulomb blockade (CB) and such system is 
named the single-electron transistor (SET). CB has been 
attracted much attention in the last decade [1-3] and 
the nature of the transport properties of SET has been 
clarified. SET is interesting because it is regarded as 
one of the most simple examples of strongly correlated 
system which can be brought into the non-equilibrium 
state by applied bias voltage. Early investigations con- 
sider the case where the tunneling conductance is so 
small that the higher order quantum fluctuation of is- 
land charge is negligible. Recently, the quantum fluc- 
tuation in SET has attracted much attention as one of 
the basic problems in this field. The quantum fluctu- 
ation is quantitatively characterized by the dimension- 
less parallel conductance: ao = i?K/((27r)^i?T) where 
^^K = h/e^ is the quantum resistance and i?T is the 
parallel tunneling resistance of the source and the drain 
junctions. There has been much development on the- 
oretical investigation in the whole range of ao- Espe- 
cially in the weak tunneling regime (ao < 1), the life- 
time broadening of a charge state level is much smaller 
than the typical level spacing of charge states, and thus 
the effective two-state model, which is equivalent to the 
multichannel anisotropic Kondo model in the equilibrium 
state[4], well describes the low-energy physics. With this 



model, it is predicted that the quantum fluctuation of 
charge causes the renormalization of the conductance 
and the charging energy below the "Kondo temperature" 
Tk = ^Q-'^/C^o'o) [4_9], The renormalization of the con- 
ductance is confirmed experimentally as 1/ In T depen- 
dence of the conductance peak at low temperature[10]. 
It is also predicted that in the non-equilibrium state, the 
dissipative charge fluctuation causes the life-time broad- 
ening of a charge state level and smears the structures of 
I-V characteristic[3, 8, 11]. 

Though investigations have revealed much about the 
quantum fluctuation, most of them have been limited to 
averaged quantities. In order to understand the nature 
of the quantum fluctuation, investigations on the higher- 
order correlation function of fluctuation operators are re- 
quired. A good starting point may be the investigation 
on the second moment of fiuctuation operators, i.e., the 
noise [12]. The charge noise and the current noise in the 
weak tunneling regime is also important for practical ap- 
plications, because it determines the performance of SET 
electrometers [13-16]. 

The current noise is defined by the auto-correlation 
function of the current fluctuation operator 51 [t) = I{t) — 

m) as 



Su{t,t') - {{5i{t),5i{t')}), 



(1) 



'Present address Max-Planck-Institut fiir Mikrostrukturphysik 
Weinberg 2, D-06120 Halle, Germany 



where (• • • ) means the statistical average. Until re- 
cently, investigations on the noise have been done us- 
ing the framework of the "orthodox" theory[13, 15, 17- 
19], which takes account of the lowest order quantum 
fiuctuation, namely the sequential tunneling (ST) pro- 
cess. Recently, several authors[14, 20, 21] discussed the 
higher order quantum fiuctuation in CB regime within 
the "co-tunneling theory" [22]. However, there is no ap- 
proximation covering both ST and CB regime. The aim 
of the present work is to construct a theoretical frame- 
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work, which covers both of these regimes for arbitrary 
ao and clarify how the quantum fluctuation affects the 



noise. 



The Keldysh formaHsm[23-25] has been one of the 

most powcrhil methods to study the non-cquihbrium 
properties of mcsoscopic systems. However to apply this 
method to SET in the two-state limit, one must overc;ome 
a technical difSculty: The spin- 1/2 operator, which is in- 
troduced to restrict charge number states by the strong 
Coulomb interaction, prevents one from utilizing Wick's 
theorem. The most successful treatment to overcome this 
problem is given in Ref. [6], in which a formulation of 
pcrturbativc expansion for the reduced density matrix in 
the real time domain is developed and the inelastic reso- 
nant tunneling processes is treated. The method of Ref. 
[6] enables one to classify various tunneling processes us- 
ing diagrammatic techniques, and can be also applied to 
other systems with the strong local correlation, such as 
quantum dot [26]. In spite of these successes, it seems to 
be still difficult to apply this method for the calculation 
of higher order correlation functions, since this method 
requires to solve a special integro-differential equation 
even for the calculation of the average in the presence of 
large quantum fluctuation [6] . 



In this paper, wc investigate the current noise and 
the charge noise in the regime of large quantum fluc- 
tuation of charge out of equilibrium. We adopt the 
Schwinger-Keldysh approach and the drone-fermion rep- 
resentation of the effective spin-1/2 operator[27, 28]. 
Schwinger-Keldysh approach enables us to calculate any 
order moment systematically by the functional derivative 
technique[29, 30] satisfying the charge conservation [29], 
and it helps us in manipulating many complicated terms. 
The drone-fermion representation allows us to utilize the 
fermionic Wick's theorem and to take effects of the strong 
correlation into account. With the help of this tech- 
nique, we can extensively take account of the higher order 
processes of tunneling. We will show that our approx- 
imation reproduces the resonant tunneling approxima- 
tion (RTA) [6] as for the average current and the average 
charge. 



The outline of this paper is as follows. In Sec. II, we 
briefly summarize the Keldysh formalism and introduce 
an approximate generating functional. We also show that 
the average and the noise expressions can be derived us- 
ing the functional derivation. In Sec. Ill, we actually 
calculate the average current, the average charge, the cur- 
rent noise and the charge noise. In Sec. IV wc show nu- 
merical results for the noise and give some discussions on 
the non-equilibrium fluctuation, the thermal fluctuation 
and the renormalization effect. Section V summarizes 
our results. 



II. KELDYSH FORMALISM AND 
GENERATING FUNCTIONAL 

A. Brief introduction of the Keldysh formaHsm 

In this section, we give preliminary definitions of the 

Schwinger-Keldysh approach and we summarize three 
useful representations: the closed time-path, the single 
time and the physical representations (Sec. 2 of Ref. 
[25]). For simplicity, we consider the following action 
of a free fermion with a linear source term to explain the 
basic formalism: 



= / dt{a{tr{ 
Jc 



ihdt - s)a{t) + a{t)J*{t) + h.c.},{2) 



where the closed time path C consists of the forward 
branch C+, the backward branch C_ and the imaginary 
time-path Cr as shown in Fig. 1[31]. a{t) is a Grassmann 
variable satisfying the anti-periodic boundary condition 
a(— oo e C-)-) = — a(— 00 — ihp e Cr)- The complex vari- 
able J{t) is defined only on the forward and backward 
branch C+ -I- C_ . The generating functional for the con- 
nected closed time-path Green function (GF) is defined 
as. 



W 



-ihlnZ, Z = J V[a*,a]exp{-S/ih). 



GF in the closed time-path representation is obtained by 
the second derivative of the generating functional with 
respect to J{t) {t e C+ + C-): 



-G(l,2) 



SJ{l)*SJ{2) 



(3) 



Hereafter, we use arguments 1, 2 instead ofti,t2 for short. 

Though the closed time-path representation makes the 
formulation compact, in order to obtain the physical 
quantities, we sometimes need the single time represen- 
tation in which the time on C is projected onto the real 
axis. In this representation the degrees of freedoms of 
fields are doubled which we denote as, 



J{t) 



J+{t) 
J- it) 



(4) 



etc. Here J±{t) is defined on C± and t is the real time. In 
the same way as J, GF is transformed into 2x2 matrix 
in the Keldysh space: 



G(l,2) = 



G+ 
G 



(1,2) 
+ (1,2) 



Q- 



-(1,2) 
-(1,2) 



(5) 



Here, arguments ti and t2 arc the real time 
and each component is defined with the statis- 
tical average in the path integral representation 
{A) = i'V[a*,a]AeyiY>{-S/ih)/Z\j^^ as G*J(1,2) = 
{ai{l)aj(2)*)\j^Q / (ih). Diagonal components G++ and 
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G~~ are the causal and the anti-causal GF, respec- 
tively. Off-diagonal components are correlation func- 
tions, which are written in the operator representa- 
tion as G-+{1,2) = {a{l)a{2)'^) / {ih) and G+-{1,2) = 
— (a(2)'''a(l))/(i?i) [25]. Here the statistical average is 
defined as (^4) = Tr[e^''^'' A]/Zo where the Hamiltonian 
operator is ea^a and the partition function Zq is given 
as Tr[e-''^«]. 

It is known that four components of Eq. (5) are not in- 
dependent. This redundancy is removed by the Keldysh 
rotation[24, 25]. After the rotation, we obtain the so- 
called physical representation: 



J = 



Ji 



G(l,2) = 



^ ^ / J+ - J_ 
V2 W+ + ^- 
G^(l,2) 
G^(l,2) G^(l,2) 



(6) 
(7) 



GFs denoted by superscripts, A, R and K are advanced, 
retarded and Keldysh components, respectively. In the 
practical calculations, instead of Ji and J2, the center- 
of-mass coordinate Jc and the relative coordinate Ja, 



Jc{t) = 



J+{t) + J^{t) 



J2(<)/%/2 



AW = J+{t) - J -it) = Ji{t)V2 



(8) 



are used in most cases. 

There are the following relations between components 
of GF: 

G— (l,2)-f G++(l,2) = G-+(l,2) + G+-(l,2) 

= i:({a(l)t,a(2)}) = G^(l,2), (9) 
G^(l, 2) - G'^(l, 2) = G-+(l, 2) - G+-(l, 2) 

= l([a(l),at(2)]) = G^(l,2). (10) 

Here we introduce a notation G*^ whose physical meaning 
is the spectral density in the energy space. Equation (9) 
is derived from the normalization of the step function 
(see Eq. (2.67) in Ref. [25]). 

The normalization condition of the density matrix re- 
sults in an important equation Z\j^q/Zq = 1 (see Sec. 
2.4 of Ref. [25]) [32] which is equivalent to the following 
equations: 



SW 



SJcil) 



Ja=0 



= • •• = 0. 



(11) 



Ja=0 



B. Model Hamiltonian in the drone- fermion 
representation and the generating functional in 
closed time path - path integral representation 

In this section, we introduce our model Hamiltonian 
and derive the generating functional for SET, based on 
which we construct the perturbation theory. 



C+ 
> — 



+00 



co-itiP 



FIG. 1: The closed time-path going from —00 to 00 (C+), 
going back to —00 (C-), connecting the imaginary time path 
Ct and closing at t — —00 — ihp. 




FIG. 2: The equivalent circuit of a SET. 



Figure 2 shows the equivalent circuit of a SET. A 
metallic island exchanges electrons with the left (right) 
lead via a tunnel junction characterized by the tunneling 
matrix element 7l(r)- The island is coupled to leads and 
a gate via capacitors Gl, Gr and Cq. We consider in 
the weak tunneling regime ao < 1 and use the two-state 
model. In this paper, we limit ourselves to the symmet- 
ric case: Gl = Gr and Tl = Tr. The total Hamiltonian 
consists of the unperturbed part and the tunneling 
part H^. Omitting a trivial constant, the unperturbed 
part is given as 



Ho 



r— L,R,I k.7i 



Ao- 



1 



(12) 



where Orfen is the annihilation operator of an electron 
with wave vector k in the left (right) lead (r = L(R)) or 
in the island (r = I) . The subscript n numbers the trans- 
verse channels including spin degree of freedom. The 
density of states is considered as constant in each region: 
Pr(£) = — £rfc) = Pr (f = L,R, I). The second 

term is the charging energy and the effective spin-1/2 
operator a acts on the lowest two charge states. The 
energy difference between two charge states is given by 
Aq = £^c(l — SQc/e) where Qq is the gate charge. 
The tunneling Hamiltonian 



r=L,R 
k,k'n 



(13) 



describes the electron tunneling across the junctions and 
simultaneous change of the charge state of the island. 
(p{t) = eVt/h is the phase difference between the left and 
the right leads and parameters kl = — '*r = 1/2 charac- 
terizes the voltage drop between the left (right) lead and 
the island. The tunneling Hamiltonian is adiabatically 
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turned on in the remote past and off in the distant fu- 
ture. It is the widely adopted procedure, which ensures 
the time translational invariance to describe a stationary 
state. 

In order to utihze Wick's theorem for fermions, we 
employ the mapping of the effective spin-1/2 opera- 
tor onto two fermion operators c and d [27, 28, 33]: 
(7+ = c'cj), (Jz — 2c^c — 1, where cf) = + d \s a. Ma- 
jorana fermion operator {(f)^ = 1). This representation is 
called drone- fermion representation [28], because is a 
"drone" whose only job is to make spin-1/2 operators of 
different spins commute, rather than anti-commute[33]. 

Employing the Hamiltonian operator in the drone- 
fermion representation and following the standard man- 
ner to introduce a path integral [34], we obtain the gen- 
erating functional in the path integral representation: 



Z = J V[a*.f.,^,a^kn,c*,c,d*,d]exp 



ih 



(14) 



where all field variables are Grassmann variables satisfy- 
ing the anti-periodic boundary condition. The action is 
given by 



dt{c{ty{ihdt - h{t))c{t) + ihd{t)* dtd{t) 



J2 T,e"^^'-'^a,knityaik'n{t)a+{t) + h.c.}, (15) 



r=L,R 
k,k'n 



The particle-hole GF a = X]r=L r "^r, m the closed time- 
path form is written as 



ar(l,2) 



z;iiVchT25r(l,2)5i(2,l)e 



j(Vr(l)-¥'r(2)) 



(17) 

where iVch is the number of the transverse channels. The 
GF for free electron in the lead r (= L, R) and the island 
r(= I), 5r(i,i') = Y.k9rk{t,t'), is given as 



9;knit,t') = {indt-e,k)S{t,t'), 



(18) 



which satisfies the anti-periodic boundary condition: 
5i.fe„(t, -oo e C+) = -5i.fe„(i, - oo). The c-field 
and the d-ficld GFs defined by 

g-\t,t') = {ihdt-h{t))5{t,t'), (19) 
g^\tX) = ihdtS{t,t')/2, (20) 

also satisfy the anti-periodic boundary condition. 

Using Eq. (16), we construct a systematic perturba- 
tion expansion in terms of a. For example, we show the 
diagrammatic representation of the zeroth {W'^^^), the 
first {W'^^^) and the second order contribution {W'^^^) to 
the generating functional W in Fig. 3. Here, solid lines, 
dotted lines and wavy lines represent GFs for c-field, d- 
field and particle-hole, respectively. Practical forms are 
given as follows: 



|Tr [{g^ll^f] , 



(21) 
(22) 

(23) 



where, ^pT:{t) = K.-[ip{t). In Eq. (15) we introduced aux- 
iliary source fields, h{t) and ip{t), in order to calculate 
the average and the noise by the functional derivation. 
It is noticed that the degrees of freedoms are doubled, as 
shown in Eq. (4). After the derivation, these variables 
are put as h±{t) — Aq and 'f±{t) — eVt/h to be related 
with the parameters of the actual system. 

By introducing a linear source term dlJ(l)(/>(l), 
where J is a Grassmann variable, all fields can be traced 
out [35]. In the limit of large transverse-channel number, 
Z is expressed as [36] 



(2) 

where W^c_ficid term corresponding to the first di- 

agram in Fig. 3 (c), which we call the c-field correc- 
tion. The other four diagrams can be written in the 
same way. Here, the self-energy of the c-field Sc(l, 2) = 
Sr=L R ^r(l, 2) is defined as 



Z 




X exp 



^ SJ 6 J 



dld2J(l)g^(l,2)J(2) 



(16) 



where we omitted the partition function of noninteracting 
electrons. The trace Tr and the products represent the 
integration along C as follows. 



Tr [5. 



dld2g,(l,2)J(2)g42,l). 



Sr(l,2) = -z?iar(l,2)5^(2,l). 



(24) 



(a) 



-\h (b) ih 




FIG. 3: The diagrammatic representation of (a) ly'"-*, (b) 
W^^"^ and (c) W'^'. Solid lines, dotted lines and wavy lines 
represent GFs for c-field, d-field and particle-hole in the closed 
time-path representation, respectively. 

In a previous paper [36], we have shown that the first 
order contribution causes the divergence at the degener- 
acy point Ag = for average charge. In order to regular- 
ize the divergence, we proposed an approximate generat- 
ing functional obtained by summing up c-field corrections 
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{qc'^c)" to infinite order: 



(25) 



Figure 4 (a) shows the diagrammatic representation of 
W. The circle represents the self-energy of c-field and 
the thick line represents the full c-field GF defined by 
the Dyson equation 



(26) 



whose diagrammatic representation is shown in Fig. 4 
(b). In Sec. IIIB we will show that W reproduces results 
of RTA for the normal metal island. In Sec. IIIC, we 
calculate the charge noise and the current noise based on 
W. 



(a) -in Q = -ih [^Q - 
O = ^ (b) 



FIG. 4: (a) An approximate generating functional including 
infinite order c-field correction. The circle is the self-energy 
of the c-field. (b) The Dyson equation for full c-field GF in 
the closed time-path representation. 



C. Formally exact expressions for average and 
noise and the charge conservation 

In this section, we summarize expressions for the av- 
erage and the noise on the basis of functional derivative. 
Relations between physical quantities in the generating 
functional representation and those in the operator rep- 
resentation are demonstrated in Appendix A. We also 
show that the gauge invariance of generating functional 
leads to the charge conservation law. 

The exact average current expression is given by the 
functional derivative of exact generating functional W 
with respect to the phase difference[29, 30]: 



lit) 



SW 



ip^{t)=eVt/hJi^{t)=Ao 



(27) 



The center-of-mass coordinate of the phase difference is 
determined by the Josephson relation[37, 38]. The rel- 
ative coordinates are put to zero because they are the 
fictitious variables. From now on, we suppress the equa- 
tions in the subscript after the vertical bar for short. The 
average charge is calculated by the functional derivation 
in terms of the scalar potential of c-field as 



Qit) 



1 



SW 



6hAit) 



(28) 



The current noise and the charge noise defined as Eq. 
(1) are given by the second derivative with respect to ip^ 
and Ha [39]: 



2e^ 



Sii{t,t') = 

LI 



ih SipAit)SipAit') 

-ihS^W 



(A 



6hA{t)ShA{t') 



4 ' 
(A-c) 



(29) 
(30) 



where (A c) is obtained from the first term by replac- 
ing subscripts A with c. 

It is convenient to rewrite Eq. (27) in the form of 
the linear combination of the tunneling current at the 
left junction /l and that at the right junction /r as 
I{t) — X]r=L R '*i-^r(0- the Same way, Eq. (29) is 
written in terms of the correlation function of It and 1^' 
(r,r' = L,R), which we denote by S'/r/i', as Sii{t,t') = 
r'=L R '^r' "^^r ^r' (t, t'). Here, /r and Sir ir' are writ- 
ten as 



Sir lT:'it,t') ~ 



5W 



h SifrAit) 

2e2 



(31) 



ih Sifr A{t)Sipr' Ait') 



(A-c) 



(32) 



by regarding ipi^ and ipn. as formally independent vari- 
ables. 

The generating functional Eq. (14) is invariant under 
the gauge transformation, i.e., the phase transformation 
of the c-field and the change of the c-field scalar potential: 



ifirit) ifrit) + 6^jit) 

hit) hit) - hSidt^Pit)) ' 



(33) 



where Sip is defined on C++C- . The relation between the 
gauge invariance and the charge conservation in the non- 
equilibrium state has been analyzed in Ref . [29] . For our 
system, the following expressions of the current continu- 
ity and the charge conservation for correlation functions 



dtQit) = ^ m, (34) 

r=L,R 

dtdt'SQQit,t') = Slrlr'it,t'), (35) 



r,r'=L,R 

can be proved (Appendix. B). 



III. APPROXIMATE EXPRESSIONS FOR 
NOISE 

In this section, we derive approximate current noise 
and charge noise expressions, which is main purpose of 
this paper. We summarize GFs in the physical represen- 
tation in Sec. Ill A which are needed for practical cal- 
culations. In Sec. IIIB we calculate the average charge 
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and the average current and show that our approxima- 
tion completely reproduces the results of RTA. We also 
give some notes on the diagrammatic rule suitable for 
the functional derivative technique. In Sec. IIIC, we 
calculate the current noise and the charge noise based on 
the diagrammatic rule. We will check that our results 
satisfy the charge conservation law and the fluctuation- 
dissipation theorem. 



A. Fourier transformation of free Green functions 

We summarize the Fourier transformation of the re- 
tarded and the Keldysh component of GFs. Hereafter 
in this subsection, we put the auxiliary source fields as 
shown in the subscript of Eq. (27). The solutions of dif- 
ferential equations (19) and (20) imposing anti-periodic 
boundary condition are [3 5] 



5^(£) = 2/(£ + i77), g$ie)=0, 



(36) 



9. 



l/{e + ir]-Ao) 



(39) 



f(£) = -2i7rtanh(^)5(e- Ao) ' ^^^^ 

where ry is a positive infinitesimal number and the 6- 
function in the energy space is defined as 

6{s)=v/{n{s' + 7j')). (38) 

The advanced component is the complex conjugate of the 
retarded component: g^(^c){£) = fjf,(c)i^)* ■ 

The two components of particle-hole GF are given by, 

a^{e) = —iT:a^p{£ — /ir) 

af-(e) = — 2i7ra°/9(e — /Xr) coth (^^^) 

where /ir = KrcV (Appendix. C). a|? is the dimension- 
less conductance for tunnel junction r written in terms 
of the tunnel resistance i?r as q° = i?K/((27r)^i?r) = 
NchT^piPr- The spectral density of the particle-hole 
propagator is given by p(e) = eE^/{e'^ + E^) where the 
Lorentzian cut-off function is introduced[6]. The Keldysh 
component and the component defined by Eq. (10) for 
the c-field self-energy Eq. (24) and those for the particle- 
hole GF are related to each other: 

Ef(£) = af(e), I]f(e) = ^.f(£). (40) 

The retarded component of the self-energy is given as 



^' 2ttT J ^ \2ttT 



2itT 

af (e) 



(41) 



where tp is the digamma function. The full c-field GF 
is obtained by solving the Dyson equation in the closed 
time-path representation Eq. (26) [36] : 

■ Gf(£) = l/(e + zr?-Ao-Ef(£)) 

Gf(£) = Gf(£){Sf(£)-2ir?tanh(^)}G^(e) ' 

(42) 



where we used the definition of (5- function Eq. (38) . Here 
we remark the following: Equation (42) shows that in 
the limit of 77 ^ 0, the charge states are independent of 
the initial equilibrium distribution, because the Keldysh 
component of c-field GF represents the distribution of 
charge states. This fact suggests that W describes a 
physically reasonable non-equilibrium stationary state, 
which should not depend on any initial state. 

We transform above expressions into single-time rep- 
resentation. Employing Eqs. (9) and (10) we obtain 



E±T(e) 



= Ta?^(£), 



(43) 
(44) 
(45) 



Here/j.^(e) — /^(e— /ir) andn~{e) = n^{e—pr) are writ- 
ten with the Fermi function f~{e) = l/(e''^ -I- 1) and the 
Bose distribution function n~(e) = l/(e''^ — !)• Func- 
tions /+ and n+ are given by /+(e) = f~{e)e^'^ and 
n+(£) = n~{e) e^^, respectively. 



B. Reformulation of the resonant tunneling 
appr oximat ion 

The reason why we adopt the generating functional 
approach is that once an approximate generating func- 
tional is obtained, one can calculate any order moment 
systematically by the functional derivation. In the follow- 
ing sections we perform practical calculations, employing 
W introduced in Sec. II B. Detailed discussions on our 
approximation are retained in Appendix D. 

The average current is calculated by substituting W 
into Eq. (31) [36]: 



-e Tr 



(46) 



Here, we used the fact that the self-energy Eq. (24) in- 
cludes the phase factor through the particle-hole GF Eq. 
(17). In the language of Feynmann diagrams, Eq. (46) 
can be rewritten in a compact form: 



(47) 



Solid dots with t represents 5(pr/5(prA{t)- The circle with 

r is the partial self-energy defined by Eq. (24). Here 
we obtain the diagrammatic rule similar to that of Refs. 
[40, 41]. 

(i) The diagrams corresponding to the functional 

derivative with respect to (p-[A(t) is obtained by a 
series of the following operations; to put a solid dot 
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onto all possible positions of vertices in a closed 
diagram, to assign r on a circle connected to the 
solid dot, to multiple i, and to assign minus sign if 
a solid line comes into the solid dot. 

Next we project the fictitious time on C to the real 
axis. As the tunneling Hamiltonian is zero on Cr, Sr(t, t') 
is zero for t E Cr or t' E Ct. Hence, the diagrams in Eq. 
(47) are rewritten as. 



Tr 



SifrAit) 



(48) 



where we performed the Keldysh rotation, (s = 
0,1,2,3) is the Pauli matrix in the Keldysh space [24]. 
The trace is carried out over the Keldysh space and prod- 
ucts represents the integration along the real-time as 



Tr [g^<fg4,] 



dld2Tr [3^(1,2)^(2)3^2,1)] 



The phase in the physical representation (fi,- is written as 
>fi{t) = iprc{t)T^ + ipr A{t)T° /2, which Icads to a useful 
relation for the functional derivative technique: 



5^,{t')/S^,A{t)^T"S{t-t')/2. 



(49) 



By using the property of the GF in the physical repre- 
sentation G{t,t'y — —T^G(t',t)T^, and that of a Pauh 
matrix t^t^t^ — — t^, we can see that the second term 
of Eq. (48) is minus the complex conjugate of the first 
term. To generalize this property, we obtain a helpful 
rule to reduce the number of diagrams: A diagram which 
is complex conjugate of a certain diagram is obtained by 
changing the direction of all lines and putting minus sign 
if the diagram includes odd numbers of vertices without 
solid dot. 

By performing the Fourier transformation, we rewrite 
Eq. (47) as: 



2eRe 



(5> 



= 2eRc 



G,(e)^S,(e)ri 



deS,^(£)Gf (e) 



-(C^K). 



(50) 
(51) 



The second term (G <-!■ K) is obtained from the first term 
by swapping superscripts K and G. By generalizing Eq. 
(50), we obtain the rule to calculate the diagram: 

(ii) Put 2x2 matrix GFs in the physical repre- 
sentation to the corresponding lines and circles, 
put T°/2 or to a vertex with or without sohd 
dot. Subsequently, carry out the trace over Keldysh 
space and the integration over the frequency e/h. 

Employing Eqs. (9) and (10), Eq. (51) 
is transformed into the single time representation: 
(e/h) J deJ:+-{e)G-+{e)-{+ ^ -). By using Eqs. (43), 



(44) and (45), and noting I{t) = h{t) = -Init), we ob- 
tain the final form of the current expression which has 
the same form as the Landauer formula: 



I(t) = 



1 



d£T^(£){/f(e)-/R-(e)} 



(52) 



where, is the effective transmission probability of lead 
electrons thorough the island: 



T^(£) 



-af(£)a^(£)|Gf(e)|2 



(53) 
(54) 



From the second form, we can see our result is equivalent 
to that of RTA[6, 42]. 

In the same way as the average current, the average 
charge is evaluated by substitute W into Eq. (28) [43], 



Qit) 



where the solid dot with t corresponds to 6h/6hA{t) 
whose practical expression is equal to the right-hand side 
of Eq. (49) in the physical representation. The diagram- 
matic rule of functional derivation in terms of the scalar 
potential is given as follows: 

(i') The diagrams corresponding to the functional 
derivative with respect to h(t) is obtained by insert- 
ing a solid dot into all possible positions of c-field 
GF. 

Following the rule (ii), the diagram for average charge 
can be also calculated: 



Qit) 



-Tr 



G.ief-^ 



de 

^^lmGf(s) 



(55) 



The imaginary part of G^, which has a peak at e ~ 
Ao, describes the excitation property of the charge state. 
When the broadening of the peak is sufficiently small and 
e, T, eV ^ Ec, G^ is approximately given by 

Gf (e) - z/(e- zAo + izImSf (zAo)). (56) 



Here z is 



the renormalization factor 1/(1 — 
a,ReEf(£)|^^^^^) ~ 1/(1 + 2aoHEc/ec)). The 
low energy cut-off ec is max(|zAo|, 27rr, |ey|/2) where 
one of three parameters must be much larger than the 
other two [11]. The imaginary part of self-energy ImE^ 
represents the life-time broadening effect. When the 
renormalization effect is negligible z ~ 1, it is written as 
ImEf (z Ao) - 7 = ftr/2, where P = Er=L,R(rir + r^i) 
and 

r.i = ^i^i^ n-{Ao - M.), Pi. = r.ie-^(^o--'). 



(57) 
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Tri (rir) is equal to the tunneling rate into (out of) 
the island through the junction r estimated by Fermi's 
golden rule. It is noticed that the condition zap <C 1 
is enough to neglect the broadening of the peak. By 
using the approximation Eq. (56), Eq. (55) for equi- 
librium state reproduces the result of RTA[6], Q/e ^ 
(l-ztanh(zAo/(2r)))/2. 

According to Ref. [6], Eq. (52) is consistent with the 
co-tunneling theory and orthodox theory for two-state 
limit[3]. We can also confirm that Eq. (55) is consistent 
with the orthodox theory in the following way: In the 
limit of ao ^ with keeping eV or T finite, the renormal- 
ization factor approaches unity and the imaginary part 
of Eq. (56) reduces to the (5-function (— I/tt) 5{e — Aq). 
Thus, Eq. (55) reproduces the result of the orthodox 
theory 

hm g(t)/e = r+/r, (58) 

where r+/_ = Er=L,RTrI/Ir- 



C. The current noise and the charge noise based on 
the reformulated resonant tunneling approximation 

Here, we calculate the current noise and the charge 
noise based on the reformulated RTA. In the following 
discussions we limit ourselves to zero frequency compo- 
nent 

5/r7r'(0) - j dtSni,-{t,t'), (59) 

where we used a fact that at a stationary state, the time 
translational invariance is satisfied and the correlation 
function depends only on the difference of t and t' . From 
the definition Eq. (29), one obtain the current noise by 
applying the rule (i) twice. Reducing the number of dia- 
grams and using the rule (ii), which can be also applied 
to the calculation of zero frequency noise diagram (Ap- 
pendix. E), we obtain 



'(0) 



2 
(-" 



2e2 
de 



di2Re 




(A-c) 



Tr 



^S,.(£)^G,(e)<5,,., + ^tAe)T'G,ie)T'tAe)YGc{e) - ^t,{e)T'G,{e)^tAs)T'G,ie) 



where we used a fact that the second term of Eq. (29) is obtained from the first term by changing t° to in the 
physical representation (Appendix. E). After some straightforward calculations, we obtain the following expression: 



5/7(0) = 5/,./r(0) --5/,.7,(0) 



G'^ie){f^is)f^is)+f^{e)f^{e)} 



G^(£) {/L(£)-/R(£)r 



2 



d£[r^ (£){/i:(e)/+(e) + /+(£)/r (e)} - {eY{f^{e) /r (e)} 



(60) 
(61) 



where r is the other side of r (for example, when the index r is L the index of the other side r is R) . The charge noise 
is evaluated by adopting the definition Eq. (30) and by applying rules (i') and (ii): 



5qq(0) = e'h' / At 



(A ^ c) _ ^2^2 f de 



4 



-Tr 



_o _o 

Ge(£)— Ge(£) — 



de|Gf(£)| af(e)a^(£)/+(£)/,7(£), 



(62) 



r,r'=L,R 



where we used Eqs. (43) and (45). 

Equation (61) has the same form as the current 
noise expression of a point contact without Coulomb 
interaction[44]. This result is anticipated, because the 



tunneling current is expressed in the same form as the 
Landauer formula. However there is an important differ- 
ence as mentioned in Ref. [11] : The effective transmission 
probability includes the Coulomb correlation and the in- 
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elastic relaxation effect. 

Here we note the following points. First, our approxi- 
mation satisfies the charge conservation law: As the ap- 
proximate generating functional W is invariant under the 
gauge transformation Eq. (33), we can show the con- 
servation law, Eqs. (34) and (35), for the approximate 
expressions. Especially for zero frequency component, 
they reduce to equations for the current conservation law, 

Sr^L.R^r = ^^d I]i.,i./=L,R '^'/r /r' (0) = 0. It is WOrth 

noticing that the gauge invariance is automatically sat- 
isfied when W consists of closed diagrams. Secondly, 
we can show that our result satisfies the fluctuation- 
dissipation theorem: At = 0, the current noise expres- 
sion is reduced to the Johnson- Nyquist formula S'//(0) = 
4TG', where G is the conductance expressed as G = 
limv^odl/dV = R^'^ J deT^ {e) / {AT cosh {e/{2T)f}. 

We discuss the current noise expression in the regime 
where zao -C 1 and the renormalization effect is negligi- 
ble, 2; ~ 1. Using Eqs. (53) and (56), the first term of 
Eq. (60) is approximately given by 



2e (7r+(7)+7r"(7)), 



(63) 



where 



7+(7) 



i?K 



de 



p{e - /(,);/, (c)/){c - /(,);), (e) 



|£-hi7- Ao|2 



(64) 



77(7) = 7+(7)e-'^'=^. 

For 7 = 0, 7i^^ (0) is the co-tunneling current from lead 
r/f to lead f/r, and Eq. (63) reproduces the co-tunneling 
theory [14, 20, 21] in the two-state limit. Equation (63), 
which we call the co-f,unneling theory with life-time broad- 
ening, is equivalent to the previously proposed equation 
in Ref. [13] in the limit of T ^ 0. The validity of this 
approximation is discussed in the next section from the 
point of view of the numerical results. 

By taking the limit ag 0, with paying attention to 
the relation lim^^o(27r7){Im[l/(e - Aq -h = ^{^ - 

Ao), we obtain the result of the orthodox theory in the 
two-state limit: 



Jim [a^'SuiO)] = a,'2e (l+ - ^ 



(65) 



where I± = e(rLiriR ± riLrR i)/r. /_ is the average 
current of the orthodox theory. The second term is re- 
lated to the part of Eq. (60) proportional to [T^Y, which 
represents the Coulomb correlation and reduces the cur- 
rent noise from the Poissonian value. It is noticed that 
even at small tunneling conductance, the second term is 
important in ST regime. 

In the same way as above discussions, for za^ <C 1 
and z ^ 1, we obtain the co-tunneling theory [14] with 
life-time broadening from Eq. (62): 



5qq(0) 



E 



e4i?K [ ^^a^{e)a^{e)f+{e)f-{e) 



=L,R 



27r2 



{(£-Ao)2+72}2 



In the limit of ao ^ 0, we can confirm that our result 
reproduces the orthodox theory [45]: 



lim [aoSQQ(O)] = ao4e2r+r_/r3 



(67) 



(66) 



IV. RESULTS AND DISCUSSIONS 

In this section we present results obtained by the nu- 
merical calculation of Eqs. (60) and (62) [46], based on 
which we discuss the non-cquilibrimn fluctuation in Sec. 
IV A and that of thermal fluctuation in Sec. IV B. We 
also discuss the effect of renormalization on the noise at 
low temperature caused by quantum fluctuation in Sec. 
IV C. 



A. Noise in the non-equilibrium state 

In order to clarify the nature of non-equilibrium cur- 
rent fluctuation, we consider the case of zero tempera- 
ture, where there is no thermal fluctuation. Furthermore, 
we limit ourselves to the condition of the high bias volt- 
age, Ec ^ |ey| ^ Tk, where the renormalization effect 
is negligible. Figure 5 shows the current noise (the top 
panel) and the charge noise (the bottom panel) as a func- 
tion of the excitation energy Aq for small ao ((a-1) and 
(b-1)) and those for large ao ((a-2) and (b-2)). Plots are 
normalized by values of the orthodox theory at Aq = 0. 
Solid lines, dashed lines and dotted lines show results of 
our approximation (Eqs. (60) and (62)), the orthodox 
theory (Eqs. (65) and (67)) and the co-tunneling theory 
(Eqs. (63) and (66) with 7 = 0), respectively. 

When chq is small, our results well reproduce the or- 
thodox theory ((a-1) and (b-1)). For large ao ((a-2) and 
(b-2)) and in CB regime (|Ao/(ey)| > 0.5), they agree 
well with the co-tunneling theory. Figure 6 shows the 
average current (a) and the average charge (b) estimated 
by our approximation (solid lines) and the orthodox the- 
ory (dashed lines). In CB regime, we can see both the 
average value and the noise are enhanced by the quan- 
tum fluctuation. Around Ao = 0, the average current 
and the current noise arc strongly suppressed due to the 
life-time broadening effect [6]. 

Next, we discuss the validity of the co-tunneling theory 
with life-time broadening (see Eqs. (63) and (66); dot- 
dashed lines in Fig. 5). As for the charge noise in the 
limit of ao — > 0, it reproduces the result of the orthodox 
theory as well as our approximation (Fig. 5 (b-1): The 
solid line and the dot-dashed line almost overlap each 
other). In ST regime it overestimates the current noise 
(Fig. 5 (a-1)), because it does not take the Coulomb cor- 
relation effect into account as mentioned before. We want 
to stress again that our result reproduces the orthodox 
theory in the limit of ao — > 0, which can not be achieved 
by the co-tunneling theory with life-time broadening. 

The physical picture of the non-equilibrium current 
fluctuation is understood more clearly with the help of 
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Ao/eV 

FIG. 5: The excitation energy dependence of the current noise 
(the top panel) and the charge noise (the bottom panel) for 
ao = 10""* ((a-1) and (b-1)) and 0.1 ((a-2) and (b-2)) at 
OK and eV/ Ec = 0.4. Plots are normalized by the value pre- 
dicted by the orthodox theory at Ao = 0: el- max = eGoV/2 
for the current noise and e^/ (4J_ max) for the charge noise, 
where Go = l/(iiL + Rn) is the series junction conduc- 
tance. The solid, dashed, dotted and dot-dashed lines show 
the results evaluated by our approximation, orthodox the- 
ory, co-tunneling theory and co-tunneling theory with life- 
time broadening, respectively. In the panel (b-1), the solid 
line and dot-dashed line almost overlap each other. The pa- 
rameters satisfy eV 2> Tk; For example, T-k/Eq ~ 10"'' for 
ao = 0.1. 



the Fano factor defined by Sn / (2e/) . The Fano factor is 
unity when the tunneUng event of electrons is Poissonian 
process and is suppressed when tunneUng events are cor- 
related. The orthodox theory predicts the sub-Poissonian 
behavior in ST regime because of the Coulomb cor- 
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FIG. 6: The excitation energy dependence of the normalized 
average current (a) and the average charge (b) for qq = 0.1 
at OK and eV/Ec = 0.4. The solid and dashed lines show the 
results of ours and those of orthodox theory, respectively. 



relation: Suppose one electron tunnels into the island 
through one junction, the next tunneling event must be 
the out going process of another electron through the 
other junction. In CB regime, the co-tunneling the- 
ory predict Sii/{2el) — 1. It means that co-tunnehng 
events, viz. the simultaneous tunneling events of two 
electrons through the two junctions, occur randomly. 

The figure 7 shows the excitation energy dependence 
of the Fano factor obtained by our approximation. In the 
small ao limit, our approximation reproduces the ortho- 
dox theory in ST regime and the co-tunneling theory in 
CB regime and smoothly interpolate two theories (The 
dotted line actually almost coincides with the result of 
orthodox theory in ST regime). For larger ao the Fano 
factor is further suppressed (the dashed and the solid 
line). Especially, our result predicts the value smaller 
than 1/2 at the degeneracy point. It is a distinctive re- 
sult because the orthodox theory predicts the inequality 
Sii/{2el) > 1/2 (p. 137 in Ref. [12]). 

Next we consider the physical meaning of our result. 
In CB regime and near the threshold voltage, the ori- 
gin of the suppression of Fano factor is considered to be 
the enhancement of the effective transmission probabil- 
ity , because the increase in the current noise is much 
larger then that in the charge noise (see Figs. 5 (a-2) 
and 5 (b-2)). As the Fano factor of the shot noise is ap- 
proximately given by 1 — T^[12], the enhancement of the 
transmission probability results in the suppression of the 
Fano factor. 
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FIG. 7: The excitation energy dependence of the Fano factor 
at eV/Ec = 0.4 for ao = 10"^ (dotted line), 0.05 (dashed 
hne) and 0.1 (sohd line). 



Around the degeneracy point, the origin of the suppres- 
sion of Fano factor is different, because the normalized 
current and the current noise are suppressed as shown 
in Fig. 6 (a) and Fig. 5 (a-2). We consider that the 
Fano factor is suppressed because of the dissipation, i.e. 
the life-time broadening effect: RTA takes account of the 
dissipation process which is the leak of an electron from 
the island while another electron tunnels into the island 
and relaxes to the local equilibrium state of the island. 
The suppression of the Fano factor by the dissipation was 
previously predicted for the ID electron channel coupled 
with a boson bath[47, 48]. 



B. Effect of thermal fluctuation 
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FIG. 8: The temperature dependence of the normalized cur- 
rent noise (a) and the normalized average current (b) at a 
threshold eV/2 = Ao = 0.2Ec for ao = lO"'^ (dotted line), 
0.05 (dashed line) and 0.1 (solid line), (c) The temperature 
dependence of the Fano factor for various qo at the threshold. 



Next we discuss the effect of the thermal fluctuation. 
Figure 8 shows the temperature dependence of the cur- 
rent noise (a), the average current (b) and the Fano factor 
(c) at a threshold for various ao- They are normalized 
by the value of the orthodox theory at T = Aq = 0. 

As the temperature increases, the average current and 
the current noise increase because of the thermal fluc- 
tuation. At sufficiently high temperature Eq ^ T ^ 
\eV\/2, orthodox theory predicts that the average cur- 
rent saturates at /_ max and the thermal fluctuation dom- 
inates the current noise, i.e. Su ~ 4(G'o/2)T which is the 
similar form as the Johnson-Nyquist noise for the ohmic 
resistance [37] (the plot for ao = 10~^ almost coincides 
with the orthodox theory). Our result further shows that 
the average current and the current noise are suppressed 
as ao increases (panels (a) and (b)). It is considered to 
be attributed to the higher order tunneling effect: The 
lifetime broadening caused by the thermal fluctuation is 
enhanced for the large tunnel conductance. The panel 
(c) is the Fano factor versus temperature plot. The Fano 
factor is independent of ao, which means that the corre- 
lation between tunneling events does not depend on ao. 



C. Renormalization effect 

Next we consider the renormalization effect at low bias 
voltage and temperature: eV,T < Tk- Figures 9 (a) and 
(b) show the charge noise normalized by c'^Rt/Ec at 
Ao = as a function of the temperature and the bias 
voltage, respectively. The charge noise is suppressed for 
large ao, which is attributed to the renormalization of the 
system parameters. In the regime eV, T <C Tk, where the 
life-time broadening effect is negligible {zag <C 1), we can 
approximate Eq. (62) as 

4(ze)2f+f_/f3, (68) 

instead of Eq. (67). Here F and r± are tunneling rate 
Eq. (57) written by using renormalized parameters such 
as Fi-i = z p{z Aq — fi,-)n^ {z Aq — /ir)/(e^i?r) where the 
renormalization factor z is 1/(1 -|- 2aoln{Ec /ec))- The 
lower cut-off energy ec is 27rT for the panel (a) and 
\eV\/2 for the panel (b). It is natural to interpret Eq. 
(68) that the charge of a carrier is modified as ze by 
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the rcnornialization effect. The interpretation is simi- 
lar to that of the doubhng of shot noise at the normal- 
metal(N)-superconductor(S) interface. Since at NS inter- 
face the carrier is 2e-charged particle, viz. a Cooper pair, 
the shot noise is twice as large as that at NN interface 
[49]. 
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FIG. 9: (a)The normalized charge noise as a function of tem- 
perature at Ao = eV — for ao = 0.1 (solid line), 0.05 
(dashed line) and 10"* (dotted line). (b)The normalized 
charge noise as a function of bias voltage at Ao = T — 0. 
Inset: The average charge as a function of Ao at T = eV = 0. 

Though the normalized charge noise is suppressed with 
increasing ao, one see that the charge noise always di- 
verges at Ao = T = = for arbitrary ao in the weak 
tunneling regime. Since the charge noise is related to 
the "charge susceptibility" for excitation energy Aq, the 
divergence means that the number of charge changes by 
"one" at the degeneracy point when we sweep the excita- 
tion energy. It is confirmed by the fact that the slope of 
excitation energy dependence of the average charge (an 
inset of Fig. 9 (b)) diverges at the degeneracy point. 

Next wc discuss the renormalization effect on the Fano 
factor. Figure 10 shows the excitation energy dependence 
of the Fano factor for various bias voltage. We can see 
at small bias voltage where the charging energy renor- 
malization is pronounced \eV\ < Tk ~ 10""^, the valley 
structures of curves are widened. The same behavior can 
be seen in the differential conductance shown in Refs. 
[3, 8, 11]. We also see that the Fano factor is suppressed 
with increasing bias voltage at Ao = 0. This suppres- 
sion is also attributed to the dissipation as discussed in 
Sec. IV A: As the bias voltage increases, the dissipative 
charge fluctuation is enhanced and thus the Fano factor 
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FIG. 10: The excitation energy dependence of the Fano factor 
for ao = 0.1 at eV/Ec = 10"* (solid line), 10"^ (dashed line), 
10"^ (dot-dashed Une) and 0.4 (dotted line). 



V. SUMMARY 

By using the drone-fcrmion representation and the 
Schwinger-Keldysh approach, we have calculated the cur- 
rent noise and the charge noise in the regime of large 
quantum fluctuations of charge out of equilibrium. We 
have reformulated and extended RTA in a charge con- 
serving way. Our approximation interpolates previous 
theories, the orthodox theory and the co-tunneling the- 
ory: Our result coincides with the orthodox theory in the 
limit of ao ^ and is consistent with the co-tunneling 
theory in CB regime. The approximation is verified 
from the fact that the result satisfies the fluctuation- 
dissipation theorem. In previous papers, we also checked 
numerically that the energy sensitivity does not exceed 
the quantum limit [50, 51]. 

We showed that at zero temperature and Ec ^ 
\eV\ ^ Tk, the life-time broadening caused by non- 
equilibrium dissipative charge fluctuation suppresses the 
current noise in ST regime. It also suppresses the Fano 
factor more than the Coulomb correlation does. Es- 
pecially the Fano factor is suppressed below the mini- 
mum value predicted by the orthodox theory 1/2 around 
Ao = 0. The origin of the suppression is attributed to 
the charge fluctuation which appears as the enhancement 
of the transmission probability in CB regime and the dis- 
sipation in ST regime. 

At Ec :$> T :$> |ey|/2 » Tk, we showed that the 
average current and the current noise deviate from the 
predictions of the orthodox theory with increasing ao- 
However, the Fano factor is independent of ao and is 
proportional to the temperature. It means that the cur- 
rent noise is dominated by the thermal fluctuation and 
the correlation between the tunneling events does not de- 
pends on ao- 

At small bias voltage and temperature eV, T < Tk, the 
charge noise is suppressed as compared with the predic- 
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tion of the orthodox theory. We showed that it can be 
interpreted as the renormahzation for the unit of island 
charge. Although the charge is rcnornialized, the charge 
noise diverges at Aq = T = eV = for arbitrary ao in 
the weak tunneling regime. It means that the quantum 
fluctuation does not wash out the charge quantization. 

In this paper, we have hmited ourselves to the discus- 
sions on the second moment and the zero frequency com- 
ponent, because wc think them primitive. Tlic investiga- 
tion on the frequency dependence of noise will be impor- 
tant to estimate the performance of high-speed SET elec- 
trometer complctely[15]. The investigation on the higher 
order moment and the full counting statistics[52, 53] will 
help us to understand carriers of strongly correlated sys- 
tem out of equilibrium. 
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APPENDIX A: RELATION BETWEEN 

GENERATING FUNCTIONAL 
REPRESENTATION AND OPERATOR 
REPRESENTATION 



which is not well defined at t = t' . Usually, an addi- 
tional operation to determine the order of field variables 
is required to remove the uncertainty. Alternatively, we 
subtract a term 



1 -ihe'^S'^W 



4 5hc{t)6hcit') 



{QAmAim 



(A3) 



from Eq. (A2) as shown in our definition Eq. (30). 
Here we used the normalization Eq. (11), {QA{t))\ = 0. 
As a result, the first term in Eq. (A2) is replaced by 
{{Q+ {t)Q- {t') +Q- {t)Q+ {t')))/2, which does not include 
the uncertainty. Using Eq. (9), we show that our defini- 
tion is equal to the standard charge noise expression 



5hA{t)ShA{t') 



(A 



{{SQit),SQit')}). 



We should stress that our definition does not change the 
final result because Eq. (A3) is from the normalization. 

The exact current expression is obtained in the same 
way. The source term corresponding to Eq. (Al) is given 

by 6S = (h/e) J dt Er=L,R^rc(i)'5¥'rA(t) + (c^ A) , 

which leads to relations 

Ir{t) = {ir{t)), STrIr'{t,t/) = {{6ir{t) , Sir' {t')}) , 

where the current operator at junction r is defined as 
Ut) = (W?i)Efc„^re'^'WaJ,„a,fc„&+ +h.c. 



APPENDIX B: CHARGE CONSERVATION 



In this Appendix, we show the relation between ex- 
pressions for average and noise in the operator represen- 
tation and those in the generating functional representa- 
tion. The variation of the exact action Eq. (15) accom- 
panied by the infinitesimal variation of h is given as the 
"twisted" combination (Sec. 9.3.2 in Ref. [25]) of Q and 
h, 



5S 



--I 



^t9MshA{t) + {c- 



A), (Al) 



where the center-of-mass coordinate of charge is Qc{t) = 
e{c+{t)*c+{t) + C-{t)*C-{t))/2. Employing this source 
term, we can show that the Eq. (28) is equivalent to 
{Q{t)) where Q = e (a^ -h l)/2 = ec^c: 



5W 



— e 



5hA{t) 



{Qcm = m)) 



where we used Eq. (9) to obtain the final form. 

The second derivative —ihe^5'^W/ 5hA{t)5hA{t')\ is 
calculated as 



{{QMQcit')) ~ {Qc{t)){Qc{t'))}\. 



(A2) 



In this Appendix, we demonstrate the charge conser- 
vation law. As W is invariant under the transformation 

Eq. (33), we obtain an identity 

-edt5W/5hA{t) = {elh) ^ 5W/5ip,A{t)- (Bl) 

r=L,R 

We can derive the other equation, which is obtained from 
above equation by replacing A with c. However, the lat- 
ter equation is not important in the following discussions. 
By putting the auxiliary source fields as the values given 
in the subscripts of Eq. (27), and employing Eqs. (27) 
and (28), we obtain the current continuity equation, Eq. 
(34). 

Next wc demonstrate the charge conservation law 
for correlation functions. By acting the operator 
iehdf S/ShA{t') or — 'ieEr'=L r ^/^'^rA(i') on Eq. (Bl), 
we obtain following two equations: 



ShA{t')6hA{t) 



9.' E 



r=L,R 



E 



ShA{t')dip,Ait)'' 

2 5^W 



r'=L,R 



SiPr'A{t') 5hA{t) ^ ^ ih Sipr'A{t') SiprA{t) 



Here the first term includes the correlation function of 4 
field variables on a same branch, {c'^{t)c±{t)c'^{t')c±{t')) , 



By comparing the left-hand side of the former equation 
and the right-hand side of the latter equation, by setting 
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the auxiliary source fields as the values given in the sub- 
scripts of Eq. (27), and by using Eqs. (29) and (30), we 
obtain the charge conservation law for correlation func- 
tions, Eq. (35). 



APPENDIX C: LOOP DIAGRAMS: 
PARTICLE-HOLE GREEN FUNCTION AND 
SELF-ENERGY 

In this Appendix, we calculate the particle-hole GF. 
We begin with the tunneling action for the large 
transverse channel obtained from Eq. (15) by trac- 
ing out the electron degrees of freedom [36]: St — 

dld2 cr+(l)a(l, 2)i7_(2). In the physical representa- 
tion, it is rewritten as / dld2 (T+(l)V^ai.(l, 2)t^(t„(2), 

where the vector field a± is defined in the same way as 
Eq. (6). Each component of can be calculated by 
utilizing the functional derivation. For example, (1,2)- 
component is 



(ri5,(l,2)ri)^ 



fc+i(l)*5a_2(2) 



-ifiTVehT^Tri 



6a* 



(CI) 



As the functional derivative is Sa±{t')/Sa±i(^2)it) = 
T°^^^ S{t — t')/^ in the physical representation, the trace 
yields e*'"'-'=^(*2-*i)/''Tr[T%.(l,2)Ti5i(2,l)]/2. Here we 
put = and ifdt) = eVt/h. The other components 



are evaluated in the same way. Then four components 
are given as 







9r 9l 

■ 9?9i 



(C2) 



where we omitted arguments and coefhcients. To obtain 
this form, we used the normalization Eq. (11) or the 
relation e{t)e{-t) = (see Eqs. (2.64) and (2.65) in 
Ref. [25]). The following calculations are same as those 
in Ref. [36]. Employing the Fourier transforms of the 
GF defined in Eq. (18), gf^ie) = -iirpr-ie) and .gf (e) = 
-2z7rpr(e)tanh(e/(2T)), we obtain Eq. (39). 

Another loop diagram, the self-energy Eq. (24), can be 
calculated in the same way. Four components are given 
in the same form as those of Eq. (C2). By using Eqs. 
(36) and (39), we obtain 



Sf(e) 



de' iaf{e') 



-, Ef(e) = a^(£). 



2Tr e + ii] — e' 

From these equations, Eqs. (40) and (41) can be derived. 

APPENDIX D: PERTURBATION THEORY 

In this Appendix, we describe some results of finite or- 
der perturbation theory and explain why we introduced 
Eq. (25) which takes account of all orders for c-field cor- 
rections. First, we calculate the first order contribution 
of average charge by employing W''^\ Eq. (28) and rules 
(i') and (ii). 




= ^(<?f(£)5f(£)Sf(£)+ffe''(^)fff(e)S?(e) + 2|gf(e)|2Ef(£)), 



where GF denoted with superscript P is given as — 
9^ + 9^ stc. In the equilibrium state, the second and 
the third terms of the second line are negligibly small 
0{Aq/ Ec) and the first term is simplified to 



tanh ( ^ ) ReEf (Ao) 



(Dl) 



where we utilized the relation g^{s)g^{£) — dAo9c{£)- 
In the limit of zero temperature, Eq. (Dl) leads to the 
log-divergence [3, 4] as -ao ln(£;c/|Ao|) sgn(Ao). 

Above result suggests that the c-field correction is re- 
sponsible for the divergence. It is further confirmed by 
calculating second order contribution of average current 



(2) 

generated from c-field correction W^c-gcid- Employing 
rules (i) and (ii), we obtain 



I^^^it) = 2eRe 




— ±,{e)T'g,{e)T'±,{e)T'g,{e) 



2eRe / ^Tr 

e / ^|gf(e)|2Ef(£)sf(e)-(C^X) + «(2) 



(D2) 



The first and the second terms of Eq. (D2), which 
are consistent of the expression for co-tunneling current 



15 



[22], also diverge at the degeneracy point. From above 
discussions, wc can deduce that the most simple way 
to regulate the divergence is to sum up c-field correc- 
tions (f^cSc)" up to infinite n as shown in Eq. (25). 
It should be noted that the correction term 51^ = 
e!Ae9^{s)Y.?{e){g'^{e)Y^^{E) - (e)}/(4/i) ~ 

0(1/77) diverges in the limit of 77 0. This divergence 
disappears when we consider Eq. (25) as discussed in 
Sec. IIIB and Sec. IIIC. 



And thus we can see that the rule (ii) can be applied for 
the calculation of the zero-frequency current noise. As 
for the second term, the derivation 



(E3) 



APPENDIX E: RULE FOR CALCULATION OF 
ZERO FREQUENCY NOISE DIAGRAMS 



In this Appendix, we demonstrate that the rule (ii) 
can be also applied to the calculation of zero-frequency 
noise. We also demonstrate the rule to calculate the sec- 
ond term of the noise expression Eq. (29) or Eq. (30). 
For example, we consider the current noise related to 
VT^-'^). Though we consider the simple case, the following 
discussions can be generalized. From the definition Eq. 
(32), the noise diagrams is obtained as: 



2e2 



= j di2Re 



(A-c) 



(El) 



The integration in terms of t in the first term is calculated 
as 



d<Tr 



Y5],(e)— 5c(e) 



(E2) 



appears instead of Eq. (49). Thus we can de- 
rive the rule that the second term is obtained 



as 



from the first term by replacing r with • 
/de2ReTr {t'^ /2)t,{e){T^ I2)gc{e) //i [54]. As a result, 
Eq. (El) is expressed as 



■I 



Re / deTr 



T°Sr(£)T05c(e) 5,,,.l{2h) - (t° ^ ri). 



As for the charge noise, we can repeat the same discus- 
sions as above. 
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